Three-Body approach to the K d Scattering Length in Particle Basis 



A. Bahaoui 1 , C. Fayard 2 , T. Mizutani 3 , and B. Saghai 4 
1) Universite Chouaib Doukkali, Faculte des Sciences, El Jadida, Morocco 
2) Institut de Physique Nucleaire de Lyon, IN2P3-CNRS, Universite Claude Bernard, 
F-69622 Villeurbanne Cedex, France 
3) Department of Physics, Virginia Polytechnic Institute and State University, 
Blacksburg, VA 24061 USA 
4) Service de Physique Nucleaire, DSM/DAPNIA, CEA/Saclay, 
F-91191 Gif-sur-Yvette, France 
(February 8, 2008) 

We report on the first calculation of the scattering length A K - d based on a relativistic three-body 
approach where the two-body input amplitudes coupled to the KN channels have been obtained 
with the chiral SU(3) constraint, but with isospin symmetry breaking effects taken into account. 
Results are compared with a recent calculation applying a similar set of two-body amplitudes, 
based on the fixed center approximation, considered as a good approximation for a loosely bound 
target, and for which we find significant deviations from the exact three-body results. Effects of the 
hyperon-nucleon interaction, and deuteron D-wave component are also evaluated. 

PACS numbers: 11.80.-m, ll.80.Gw, 13.75.-n, 13.75. Jz 
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While the threshold behavior of the A-nucleon system has been found simple, the corresponding one for the K- 
nucleon (KN) is quite complicated as its threshold is above those for the ttY (Y = A, S) channels to which it couples 
strongly Besides, it also couples to the below-threshold A(1405) resonance. Moreover, this topic has suffered 
from years of persistent ambiguity in the sign of the real part of the K~p scattering length a^-p- sign from the 
scattering data opposite to the one from kaonic hydrogen atomic data. Even under this circumstances, a few three- 
C*~) ■ body calculations on the K~ - deuteron scattering length A K - d were performed with different degrees of refinement, 
by always disregarding the controversial kaonic hydrogen constraint on Re(a^- p ) Some of these works were 

devoted primarily to the calculation of the mass and momentum distributions such as m(iTY), in the breakup reactions: 
K~d — ► ttNY , so the K~d scattering length was, to some extent, a by-product jl|j|]. Calculations required various 
two-body amplitudes as input, the most important of which being the coupled K N, ttY channels. Those amplitudes 
"--^ ■ were derived from ad-hoc rank one separable potentials with (energy independent) strengths, and ranges in the form 
»^ \ factors determined by fit to the low energy K~p scattering data. On the average the thus obtained values for A K - d 
were centered around ~ (—1.5 + il.Q) fm. Due to the very restricted quantity and quality of the data and to the 
lack of sound theoretical guidance (apart from isospin symmetry) on the form of the potentials, along with the then 
troubled Re(a,K- p ), it appeared meaningless to continue this theoretical endeavour any further. So, the investigation 
in the subject became dormant. One very important finding, however, was that the iterative solution for K~d did 
^ , diverge, hence solving the three-body equations without truncation became a must. 

Recently, there has been a steady progress in effective low energy hadronic methods such as Chiral Perturbation 
Theory flqfl. This advance as well as the new K~p data have created a renewed interest in the physics with low 



energy kaons, to the extent that there have even been discussions on extracting the kaon-nucleon a terms which are 
expected to provide important information on chiral symmetry breaking, strangeness content of the nucleon, etc. 
[Ecj-O]. Note that both a^- p and A K - d are vital ingredients in this respect. 

On the experimental side, the long-standing sign puzzle in a^- p got finally resolved by the KEK X-ray measurement 
in the kaonic hydrogen HI . The extracted scattering length is: a K - p = (-0.78±0.15±0.03)+i(-0.49±0.25±0.12) fm. 
Though the sign of the real part is now settled, one clearly needs a more accurate value, particularly for its imaginary 
part. With this in mind, remeasuring this quantity along with extracting A K - d from kaonic atom experiments is 
underway in DEAR experiment at DA<i>NE, see e.g. This should, in principle, allow for an extraction of the 

scattering length a#-„ (see e.g. Ref. Q). 

Interest in improving the calculation of A K - d may be witnessed in two recent publications. First, Deloff |l5| ] 
compared the results of old generation multi-channel three-body calculations jjj] with a simplified three-body result 
keeping only the K~p and K~n, and AA(deuteron) input (all in S- wave), and with the fixed centre approximation 
(FCA) applied to the simplified three-body model. Here the positions of the proton and neutron in the target deuteron 
were frozen at a certain separation while the KN amplitudes were replaced by their scattering lengths. The K~d 
amplitude was then obtained algebraically as a function of the proton-neutron separation. To include partially the 
effect of Fermi motion, its expectation value over the separation was calculated with the deuteron wave function, (this 
leads to the results called FCA-integ in Second, Kamalov et al. jl6| performed yet another FCA calculation, 

but with an essential difference: the input ATA potentials for the 5- wave amplitudes were obtained at 0(1/ f 2 ), lowest 
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order in the SU(3) non-linear chiral Lagrangian for the interaction of the pseudoscalar meson and l/2 + baryon octets 
|^7| . Only two free parameters were involved: the best fit to the data was found with a cut-off in the momentum 
integration at p m ax = 660 MeV, and with an effective meson decay constant / only 15% larger than the physical 
pion decay constant: = 93 MeV. With the hadron physical masses resulting from the isospin symmetry breaking, 
which the authors called the Physical (or Particle) Basis as compared with the Isospin Basis, the obtained amplitudes 
for the coupled KN, wY, rjY channels allow to reproduce the existing low energy data quite well (for References, see 
fljl). The A(1405) resonance was also generated as a bound state below the K~p and if n thresholds. This approach 
is in sharp contrast to the models mentioned earlier, in which the only constraint on the numerous parameters was 
the x 2 ^ to the available cross sections, etc. (note, however, that improved models of this type exist with SU (3) 
constraints on the relative strengths of the potentials |Q ) . 

Here we have chosen to employ a strategy similar to the one in Ref. [jlTj for determining the essential part of the 
input to the three-body equations, and solve them exactly. In this way, we will be able not only to provide the 
best theoretical value for A K - d to date, but also to test the reliability of the FCA, the effect of the irN and YN 
interactions, etc. on this quantity, as investigated in [ fL5[ within the old scheme. 

We have introduced two distinct sets of potentials which are slightly different from the one in jL7|. The main 
reasons are: (i) to check the sensitivity of the calculated A^-d to the two-body input within a reasonable margin of 
difference, and (ii) to embody them in our current investigations on the finite energy K~d scattering including the 
three particle final states like irNY , for which the momentum integration must be done along a rotated line in the 
complex plane. For this objective, instead of truncating the integration at p m ax, the potentials should have a smooth 
cut-off by form factors. Following closely Eqs. (1) to (9) of [0, the first set of potentials (= OS1) is expressed, using 
the isospin notation, as: 

V(I)a = -^Clgip^iei + ej)g{pj), 

where pi and are the magnitude of the center of mass momentum and the corresponding meson energy in the i-th 
channel, respectively. The SU(3) coupling coefficients are C^ =0 = Dij and , Cjf 1 = Fij, as defined in Tables II and 
III of |jl7). The form factor has been chosen as g(p) — f3 2 /(p 2 + (3 2 ) for all the coupled channels. A fit to the data 
with comparable quality to JljJ has been reached with /3 = 870 MeV (4.41 fm -1 ) and / = 1.20/^. The second set 
of potentials (= OS2) introduces the possible SU (3) breaking effect in the coupling strengths such that its form is 
identical to the one for OS1, except that it is now multiplied by an extra coefficient b\y By performing a standard 
statistical fit to the data, we have obtained f3 — 865 MeV (4.39 fm _1 ) and / = 1.16/ w . The values of the SU(3) 
breaking coefficients all stay within 20% around unity, see Table |. Note that, unlike in |18[ |, the radiative capture 
K~p — > jY has not been investigated. Overall, the fit to data by these two interactions and the one in |L7j are just 
about the same: differences may be examplified in terms of the scattering lengths shown in Table [l]|. All of them have 
been evaluated at the K~p threshold (=1432 MeV): beware the discussion below regarding the value of the threshold 
at which these quantities are calculated. As compared with experiment p3| , both the real and imaginary parts of 
a K-p given by all models adopted are found within 2a s t a t of the central values. The extra parameters in OS2 make 
the results somewhat distinct from the two other models. The symmetry breaking effect in the mass of the hadron 
isospin multiplets on the scattering lengths is quite visible, especially on the real parts, as one can see in Table O. (In 
the limit of isospin symmetry, one has a p = a°, and a ex = a p — a n ). Finally, we should note that just like in [ [17] we 
have retained also the i]Y channels to obtain a reasonable fit to some data like the 7rE mass spectrum. 

Other two-body input for our three-body equations consists of the TV TV interaction in the deuteron channel, and 
the P33 ttN and S'-wave YN interactions. The first one not only holds the initial- and final- state proton and neutron 
to form a deuteron, but supplies the NN scattering in the presence of a spectator kaon in intermediate processes. We 
have adopted a specific model of our choice, but taken also two other models for comparison, as discussed later. The 
remaining interactions have been taken from HQ]. 

In our three-body calculation, we first retain two-body K N ^-matrices only: for the elastic K~p, K~n, K n, and 
charge exchange K~p <-> K n, which is in line with pll. It turns out that, with only these two-body channels for 
K~d at threshold, effectively there is no other strong branch cut along the real axis in the momentum integration in 
the three-body equations, so no contour rotation into the complex plane is needed for integration, and even a sharp 
cut-off may be imposed. Thus with the two-body input fro m |17j ], we were able to find the exact solution to the 
three-body equations without making the FCA as adopted in Jl6[. Table [II summarizes our calculations. The result 



with the amplitudes from |17| is presented as in column Oset-Ramos, along with our own sets of two-body input OS1 
and OS2. For later discussions we have separated the results into: (i) pure elastic case: K~ multiple scattering on 
the proton and neutron, (ii) the total contribution, (iii) the intermediate charge exchange contribution, which is the 
difference between the values in (ii) and (i). 
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We first present the consequences resulting from the on-shell properties of the KN input. Although not entirely 
free from ambiguity, one possible way to define the on-shell contribution may be given by the corresponding FCA 
result. Given that the deuteron is very loosely bound, this could in fact be guaranteed by Beg's theorem p9| : if the 
ranges of interactions for the projectile and target constituents between two successive collisions do not overlap, the 



projectile-target interaction is described entirely by the on-shell properties of the two-body input. From Table II], 
we see that the results for all three models are more or less the reflection of the differences in the scattering lengths 
in Table jnj Now there is a bit of trouble in the present situation: near the threshold the K~p and K n elastic, and 
K~p <r+K n charge exchange amplitudes all vary rapidly due to the proximity of the A(1405) resonance. In fact, the 
minimum of the real part of the K~p amplitude is found to be located slightly below the K~p threshold (Wth = 1432 
MeV), see e.g. Fig. 9 of flijj . Besides, the threshold is slightly different for each physical K N channel, except in the 
limit of exact isospin symmetry. So, depending on the threshold energy adopted in determining the scattering lengths 
a^-p, a^° n and a K - H j» for use in the FCA, the resulting on-shell contribution to A K - d has been found to vary 
up to at least 20% for its real part. On the other hand, its imaginary part is relatively stable. It may be useful to 
remark that this strong variation in the present FCA result is due to the violation of Beg's theorem: the finite life 
time of the A(1405) causes its propagation, hence non-overlapping of the interaction ranges does not materialize. 



Now, we wish to underline a significant finding of the present work: as one can see in Table III, the Faddeev results 
for all three models are closer to each other than in FCA, with Re(A K - d ) for OS2 only about 15% different from 
the values given by the other two models. By comparing the exact three-body result and its FCA version for a given 
set of two-body KN interaction, there is a noticeable difference which may be regarded as due to off-shell effects. 
Particularly, the effect of the charge exchange scattering K~p <-> K n in multiple scattering process has been found 
grossly overestimated in the FCA. This is because the K n channel has a higher threshold than that for K~p. The 
constant scattering length approximation adopted in FCA ignores this aspect. Within the FCA, the situation gets 
even worse with the Isospin Basis in which the two thresholds are identical: see e.g. Table II of p6| . 

Next we have checked the dependence on deuteron models. First, we have compared the result with our [] SF(6.7) 
deuteron |2(J with the one using the relativized version of the model elaborated in The difference in A^-d was 
found mostly in the imaginary part, but was only within a few percent. But when a simple 3 Si-wave model is used, 
this difference grows to be about 20% as seen in the second and third columns of Table However, the real part 
appears quite stable. The short range part of the deuteron wave function should be responsible for this difference, 
and at least one needs to retain a realistic model with the 3 Di component. In the same Table, we give a comparison 
between the results in the Particle and Isospin bases. The difference is remarkable in the imaginary part, which may 
be easy to understand since all the KN thresholds are identical in the Isospin Basis, hence charge exchange scattering 
is kinematically "elastic", as discussed in the last paragraph. 



We then want to check the claim in |15| that the FCA is rather reliable relative to the full three-body result. In 
fact, by comparing the rows for FCA-integ and Faddeev in Table II of |l5|, the author seems to be right: the two 
methods provide almost identical imaginary parts, while the FCA tends to slightly underestimate the magnitude of 



the real part. This is just opposite to what has been found above: see Table HI. Eventually, we have solved this 
apparent puzzle: by taking a pure iS-wave deuteron and also by excluding the charge exchange contribution in the 
KN input to the three-body equations, we have found that the exact and FCA solutions present very similar values 
for the imaginary part, but that the latter underestimates the real part by about 30%. In fact this is how the author 
of pE| performed his calculation, and the characteristic of the outcome was just the same: main difference in the real 
part. Then, once the charge exchange contribution is introduced, we find that the trend changes considerably. We 
have found further that by introducing a realistic deuteron with the _D-component, even the result without charge 
exchange process does not satisfy the finding of |1| . Hence we conclude that the FCA is not as reliable as claimed in 

Lastly, we need to check the effects due to the irN and YN interactions, which have been excluded so far from our 
two-body input: they introduce the tt(YN) and Y(irN) states in the three-body equations, where particles outside 
the parenthesis are the spectators. Our preliminary results show effects smaller than 5%, so the semi-quantitative 
estimate of jl6| seems justified, thus the claim towards the end of jlf| appears on no solid ground at this point. 

To summarize, starting from a study of K~p scattering length, reproducing well enough the data, we presented a 
relativistic Faddeev approach for the K~d scattering length and investigated its sensitivity to various input ingredients. 
The obtained values for A K - d agree with each other within ±20%, leading to A K - d « (—1.8 + £1.5) fm. Here, our 
approach embodied elastic and inelastic KN channels in the three-body formalism. To go further, we are in the 



1 The parameters are fitted to the static properties of the deuteron, with D-state percentage value Pr> = 6.7%, and to the 
monopole charge form factor up to ~ 6 fm -1 . 
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process of including all other relevant inelastic channels, such as irY and r]Y. How one may extract the scattering 
length ax~ n from the experimental values of a^- v and A K - d , is an other question under study. A more extensive 
account will be reported in a forthcoming paper. 
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TABLE I. SU-(3)-symmetry breaking coefficients 6y (= 6^) for model OS2. 
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TABLE II. KN scattering lengths (in fm) calculated at W = M K - 


+ M p in the particle basis with models OS1 and OS2. 


The values in the last column have been evaluated by Ramos [17] at the same energy. a p , a n , a° , 


and a EX are the scattering 


lengths for elastic K p, K n, K n, and char; 


ge exchanj; 


;e if p ^K n 


, respectively. 




OS1 






OS2 


Oset-Ramos 


a p -1.04 + i 0.83 






-0.71 +i 0.92 


-1.01 + i 0.95 


a n 0.57 +i 0.45 






0.71 +i 0.69 


0.54 + i 0.53 


a° n -0.60 +i 0.89 






-0.23 + i 0.97 


-0.52 + i 1.05 


a ea: -1.37 + i 0.48 






-1.16 + z0.39 


-1.29 + i0.48 



TABLE III. K~d scattering length (in fm) calculated in the particle basis, with the FCA approximation, and with the 
Faddeev three-body model. The FCA and the Faddeev calculations in the last column have been performed by us with the 
Oset-Ramos KN scattering lengths given in Table |n[ 



FCA 


OS1 


OS2 


Oset-Ramos 


el. only 


— 1.32 + i 1.10 


-1.09 + i 1.41 


-1.36 + i 1.26 


charge ex. 


-0.83 + i 0.82 


-0.64 + 10.35 


-0.63 + i 0.69 


total 


-2.15 + i 1.92 


-1.73 + i 1.76 


-1.99 + i 1.95 


Faddeev 


el. only 


-1.70 + i 1.31 


-1.41 +i 1.48 


-1.68 + i 1.33 


charge ex. 


-0.29 +i 0.34 


-0.27 + i0.18 


-0.24 + i0.25 


total 


-1.99 + i 1.65 


-1.68 + 11.66 


-1.92 + i 1.58 


TABLE IV. 


K d scattering length (in fm) calculated with models OS1 and OS2, in the isospin basis, and 


in the physical 


basis. SF(6.7) and 3S1 specify the deuteron model. 






Model 


iso-SF(6.7) 


phys-SF(6.7) 


phys-3Sl 


OS1 


-1.76 + i 2.91 


-1.99 + i 1.65 


-1.98 + i 1.31 


OS2 


-1.37 + i 2.68 


-1.68 + i 1.66 


-1.69 + i 1.33 
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